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ABSTRACT 

A semi-empirical model is described for predicting unsteady aerodynamic forces on 
arbitrary airfoils under mildly stalled and unstalled conditions. Aerodynamic forces are modeled 
using second order ordinary differential equations for lift and moment with airfoil motion as the 
input. This model is simultaneously integrated with structural dynamics equations to determine 
flutter characteristics for a two degrees-of-ffeedom system. Results for a number of cases are 
presented to demonstrate the suitability of this model to predict flutter. Comparison is made to the 
flutter characteristics determined by a Navier-Stokes solver and also the classical incompressible 
potential flow theory. 
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NOMENCLATURE 


Aq, Aj, etc 

coefficients in the model. Equations ( 1 ) and ( 10 ) 

a, a h 

distance of elastic axis from the origin 

^oo 

speed of sound 

»1 

coefficient in the model. Equations ( 1 ) and ( 10 ) 

b 

semichord 

c h 

damping coefficient in plunging 

Cl 

lift coefficient 

C M 

moment coefficient about elastic axis 

Ca 

damping coefficient in pitching 

c 

chord 

h 

plunging (bending) displacement, positive downward 


polar mass moment of inertia 

K h 

plunging spring constant 


pitching spring constant 

M 

Mach number 

m 

mass per unit length 

Qh 

generalized force in plunging 

Qa 

generalized force in pitching 

Re 

Reynolds number based on chord 

r a 

radius of gyration about elastic axis 

S 

static unbalance, mbx a 

s 

nondimensional time, Vt/b 

t 

time 

V 

resultant air velocity 

v F 

flutter velocity 

y* 

nondimensional velocity V/bco a 

W 0.75C 

upwash at 3/4 chord 

w 

nondimensional upwash, w 0 75 C /V 

x a 

distance between elastic axis and center of gravity 

a 

pitching (torsion) displacement, positive nose up 
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a s 

<l>(s) 

M- 

Ca 

Ch 

P 

a 

x 

CO 

co A 

co h 

®a 


steady angle of attack 
Wagner function 
mass ratio 

damping associated with fluid forces 
critical damping coefficient in plunging 
critical damping coefficient in pitching 
air density 

time scale for superposition integral 

nondimensional time, t a^ c 

small amplitude harmonic oscillation frequency 

frequency associated with fluid forces 
uncoupled plunging frequency, VKj/m 

uncoupled pitching frequency, VK^/Ia 


3 



INTRODUCTION 

The aeroelastic response and flutter analysis of a helicopter rotor or an airplane wing or a 
turbomachine blade requires accurate information about aerodynamic forces under stalled and 
unstalled conditions. The theoretical approaches for prediction of aerodynamic forces including the 
effects of flow separation consist of the Navier-Stokes methods, the discrete vortex methods and 
the zonal methods. In the Navier-Stokes methods, the fundamental equations of fluid mechanics 
are solved by numerical techniques The discrete vortex methods solve for potential flow 
without the boundary layer ^>5,6 The viscous nature of the flow is modeled by the generation and 
subsequent induced transport of discrete combined vortices. The manner and location of their 
generation is obtained empirically or through appropriate boundary layer calculations. The zonal 
methods model the flow as a combination of viscous, inviscid and transition zones 7 ’^. In the 
numerical implementation of the model, these zones interact in an iterative manner. These 
theoretical models lack generality (due to various assumptions and restrictions of the formulation) 
and are computationally expensive. They are extremely difficult and expensive to incorporate into 
aeroelastic analysis codes. 

For use in a standard aeroelastic analysis, various empirical and semi-empirical models 
have been developed to model the nonlinearity and phase changes in airfoil forces with respect to 
its motion. These methods are computationally fast and allow for various airfoil motions and flow 
conditions. These phenomenological models are not based on the first principles of fluid 
mechanics, but are obtained from understanding the physical relationship between the forces on 
airfoil and its motion. The different ways for modeling the nonlinearity in lift and moment are (1) 
corrected angle of attack, where the actual angle of attack is corrected to get the effective angle of 
attack 9, which is a very popular approach in the aerospace industry, (2) time-delay, synthesis 
procedures, which involve curve-fitting experimental data with algebraic or transcendental 
functions and (3) ordinary differential equations, where the relationship between lift and angle 
of attack is described with an ordinary differential equation along with parameter identification 
from experimental data * 

Tran and Petot 1 1 originally proposed the idea of using ordinary differential equations to 
model the unsteady behavior of airfoil loads. The form of the ordinary differential equations was 
obtained by curve-fitting the experimental data for various airfoils. Substantial manual adjustment 
was needed to that model to predict accurately the airfoil loads. In the present work, a different 
semi-empirical model is developed to predict the lift and moment loads on an airfoil using an 
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ordinary differential equation. The form of equation is motivated by physical considerations. It is 
assumed that the fluid lift force has the characteristics of a damped harmonic oscillator because of 
the fluid’s inertia, dissipation and compliance. The airfoil motion is assumed to drive the fluid lift 
oscillator and the driving terms have nonlinear as well as linear components. The pitching moment 
is also represented using the same form of equation. The parameters in these equations are 
determined from comparison with airfoil load data for small amplitude oscillatory motion about a 
large steady angle of attack from a Navier-Stokes solver The model is then compared with the 
Navier-Stokes solver for large amplitude oscillatory motion of airfoil. 

This model is then used to provide the aerodynamic loads to a two degrees-of-freedom 
structural dynamics model. The coupled second order ordinary differential equations are solved 
using an implicit linear multistep method with backward differentiation formula, also known as 
Gear's stiff method. The time response for a two degrees-of-freedom structural-fluid system using 
these equations is compared with those from the Navier-Stokes solver and the classical 
incompressible potential flow model. 


MATHEMATICAL AND NUMERICAL FORMULATION 

It is postulated that the lift coefficient of the airfoil, Cl, is related to the motion of the 
airfoil, its angle of attack, a, through an ordinary differential equation of the following form: 

Cl + 2 Ca(oaCl + a>& Cl= A 0 + Ai a +A2 a 2 + A3 a 3 + Bi (a + a) (i) 


# 

where ( ) denotes a time derivative. 

This equation at present is postulated and not derived from first principles, e.g. from the 
Navier-Stokes equations. On the left hand side it is assumed that the fluid lift force has the 
characteristics of a damped harmonic oscillator because of the fluid's inertia, dissipation and 
compliance. On the right hand side, the body motion is assumed to drive the fluid lift oscillator and 
the driving terms have nonlinear as well as linear components in a, a and a . The a, a and Cl 
terms allow for modeling the phase shifts between Cl and a. Once £a> ®A’ Aj, A 3, are 

known, it is relatively simple to solve Equation (1) by analytical or numerical means for any 
arbitrary time dependent motion in combination with the structural dynamics equations. £a an d 0 >a 
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can be physically interpreted as a characteristic damping and frequency of the fluid. They are in 
general functions of a. 

To determine £ A , co A , Aj, A 2 , A 3 , and Bj for a given airfoil at some Mach number and 

Reynolds number, a numerical or physical experiment is undertaken. First in a static experiment, 
Equation (1) simplifies to 


C Ls = ( A 0 + A 1 a s + A 2 a s 2 + A 3 a s 3 ^ f W A 2 (2) 

where the subscript s refers to static or steady flow conditions. Note that co A 2 appears for static 

conditions in Equation (2) due to the choice of harmonic oscillator representation in Equation (1). 
By curve-fitting the experimental data of C Ls vs a s , the coefficients A ]/co A 2 , A 2 / 00 A 2 and 
Ayco A 2 are determined. 


Now the airfoil is oscillated about a c i.e., 
a(t) = a s + a(t) 





a 

« 

a s 


To determine £ A , co A , B), it is sufficient to consider | a | << | a s | , i.e. only (infinitesimal) small 

oscillations about various static angles of attack need be considered to determine all the model 
parameters. In analogy to Equation (3), the lift coefficient can be written as, 


Cl - Cls + Cl 


(4) 


For a « | a s | , then Cl « | ClJ and Equation (1) becomes 


Cl+ 2 Ca^aCl + w^Cl = (Ai + 2A2<x s + 3A3a£) a + B i( a + a) (5) 

Note that Equations (3), (4) and (5) should hold for both attached and separated flow for a 
sufficiently small, i.e. one can still linearize about a steady, separated flow state. This linear 
relationship significantly simplifies the determination of £ A , co A , Bj. For small amplitude 

harmonic oscillations, substituing 
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a = ae i ®t,C L = C L e i ®‘ 


( 6 ) 


in Equation (5) gives. 


C L = a 


(A1+2A2 a s +3A3 a| - Bj to 2 ) 2 + (Bi co) 2 


(0^- to 2 ) 2 + (2 C a co a co) 


m 




where 


§ = tan' 1 


Bjco 


Ai + 2A2 ots + 3A3 a| - B 1 co 2 ] 


- tan 


-1 


2£a co a ft> 


L o^-co 2 


Considering the real part of Cl Equation (6), 


(7) 


( 8 ) 


C L = 06 


(A1+2A2 (Xs+ 3A3 « 2 - Bj co 2 ) 2 + (Bj w) 2 

(co^-co 2 ) 2 +(2CaO)a®) 2 


1/2 

COS(G)t+<(>) 


(9) 


/S /N 

^A» an<1 B 1 are determined by fitting Equation (9) to Cl vs. a data from the Navier-Stokes 
solver. After determining etc, Equation (1) can be used to determine Cj^(t) from a given a(t) 
even when the motion is large and transient (not just small and harmonic in time). 


The moment coefficient of the airfoil, C^j, as a function of its angle of attack, a, is also 
modeled in a similar manner as follows 


Cm + 2CacoaCm + © i Cm = A 0 + Ai a +A2 a 2 + A3 a 3 + Bi (a + a) (10) 

The coefficients ®a> ^b ^2’ ^3» an(1 Bj for Cjyj are determined using the same 

procedure and the data from the Navier-Stokes solver. These coefficients may have different values 
for Equation(l) and Equation(lO). On physical grounds, one may expect (hope) that £a ancl W A 

will be similar for Equations (1) and (10). 
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Note that the coefficients of the semi-empirical model were determined from a combination 
of theoretical/numerical and physical/experimental data for both attached and separated flow. The 
Sankar code was used for a < 12° to obtain both dynamic and static aerodynamic data. At higher a 
the Sankar code displayed oscillating lift even for a steady angle of attack and the results appeared 
sensitive to time step, artificial viscosity, grid spacing and other parameters. Therefore to determine 
the static coefficients, e.g. Aq, Aj, etc., for 12 ° < a < 22 ° the experimental data from 

McCroskey 12 were use d. The available dynamic experimental data did not appear sufficient to 
independently assess the dynamic characteristics of the semi-empirical model and were not used. 

It should be noted that according to the Sankar code, static separation begins at about 
a = 10° and at a = 15° the separation is still mild. 

The Navier-Stokes solver was run on a CRAY XMP computer to generate small amplitude 
oscillation loops for lift and moment on a NACA 0012 airfoil at a Mach number of 0.3 and a 
Reynolds number of 4 million. The mean angle of attack was varied between 0° and 12° with the 
amplitude of oscillation being 0.1°. These values of Mach number, Reynolds number and angle of 
attack were selected, as the Navier-Stokes solver was validated over this range of values and gave 
physically meaningful results. Equation (9) was fitted for mean angles of attack between 0° and 
12°. Equation (2) was fitted to physical/experimental data for static angles of attack between 0° and 
15° to incorporate the static stall behavior into the model. If a Navier-Stokes solver were available, 
that would generate reliable data in dynamic stall regime, the proposed model could be easily 
extended to predict dynamic stall behavior with reasonable accuracy. 

The values of co^, Aj, A 2 , A 3 , and B j for and Cjyj are given in Appendix 1 for a 

NACA 0012 airfoil at M = 0.3 and Re = 4 million. These parameters were determined using a 
nonlinear regression fitting procedure from the Statistical Analysis System (SAS) A 
comparison between this model and the Navier-Stokes solver for a large amplitude oscillation is 
shown in Figures 1 and 2. It can be seen that the model, which was fitted to small amplitude 
oscillation data, predicts large amplitude oscillation data accurately. These results are typical for the 
range of angles of attack considered in this paper. 

To incorporate a plunging motion of the airfoil into the model, a is replaced by (oc + h/V) 
in Equations (1) and (10). This model is now integrated with structural dynamics equations for a 2 
degrees-of-freedom system to analyze the flutter behavior of a NACA 0012 airfoil. 
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STRUCTURAL DYNAMICS MODEL 


A typical section model of an airfoil, as shown in Figure 3, is used for the flutter study. 
The typical section has two degrees-of-freedom, plunging (h) and pitching (a), positive as shown. 
The governing equations for this model are 


mh + S a + Chh + Khh = Qh 

Sh + I a a + C a oc + K a a = Q a (11) 

where m is the mass, S is static unbalance, Cj,, C a are structural damping parameters, and K a 

are the pitching and plunging spring constants, and Q a are the generalized forces, and ( ) is 

d( )/dt. Defining h = h/2b, x a = S/mb, r& = Ia/mb 2 , a$ = Kf/m, G)a = Ka/I a , Ch = 2Ch ®h m 
and C a = 2Ca a) a I a , where b is the semi-chord, I a is the mass moment of inertia about the 
elastic axis, coj, and 0) a are the uncoupled plunging and pitching frequencies, Ch and Ca are the 
critical damping factors in bending and torsion, Equation (1 1) can be written as 

h+[^] d + [2Ch<Qh]h+[(ogh = ^- 


[M]h + [!k] 2 a + 


12 


Ca ®a rii 


a + 


rq 


a 


_ Qa 
4mb 2 


( 12 ) 


The generalized forces, and Q a are related to the lift and moment coefficients as 
follows: 


Qh=-lpV 2 (2b)CL 

Qa = 2 P V 2 (2b) 2 Cm ()3) 

where is the lift coefficient, and Cjyj is the moment coefficient about the axis of rotation (elastic 
axis), V is the resultant air velocity and p is air density. Defining, t = t aoo/2b, and V* = V/(b 0) a ), 
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where aoo is the speed of sound and substituting Equation (13) into Equation (12), the final 
aeroelastic equations can be written as 


h + 


m a+ 

4 ChMcOh 

h + 

2Mt0h 

» 

< 

* 

# 

1 

LV*gJ 

L 2 J 



12 2ClM 2 

Ttjul 




U i&M 




a + 


Mr, 


QL 


L V* 


I 2 a _ Cm M 2 
7C|i. 


(14) 


where ( ) is now d( )/dt , M = V/a^ and |X is vajn p b . Equations (13) and (14) are used with 
any one of several aerodynamic theories. 

Equations (1) and (10) for the semi-empirical aerodynamic theory are rewritten to include 
the effect of plunging displacement as follows: 

Cl + 2Ca®aCl + co& Qj= A 0 + Ai 0 +A 2 0 +A 30 +Bj (0+ 0) 


Cm + 2Ca(oaCm + Cm = Ao + A i 0 +A 2 0 2 + A 3 0 3 + Bi (0 + 0 ) 


(15) 


where 0 = a + h/V. Substituting h and t for h and t, and 0 = a + h/M. For simplicity, the bar on h 
and t will be dropped for rest of the paper. Strictly speaking a should be replaced by a + h/V 
Equation (1), but a and a should be left alone . (However, for illustration purposes the present 
approach should suffice). 

The Equations (14) and (15) form a system of four coupled, second-order ordinary 
differential equations. These equations are written as eight simultaneous first-order differential 
equations with corresponding initial conditions. These equations have widely varying time 
constants and can be described as a stiff system. The explicit methods, such as Runge-Kutta 
method, require a large number of small time steps in order to satisfy the accuracy requirement for 
the global error in a stiff system. So this initial value problem is solved using an implicit linear 
multistep method with a backward differentiation formula (Gear's method) from the IMSL 
library 
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THE NAVIER-STOKES MODEL 


The Navier-Stokes solver, that was used to generate data for fitting the semi-empirical 
model, was also used to provide aerodynamic loads for time-marching the solution to Equation 

(14). This code solves the full Navier-Stokes equations for 2-dimensional, unsteady, viscous 
flows in a body-fitted moving coordinate system. The code calculates and Cjyj for initial values 

of a and h for an airfoil. These values are then used to solve Equation (14) for new values of a 
and h at the next time step. These new values of a and h are then used by the code to calculate new 
Cl and Cjyj and so on. The fluid-structure interaction is thus analyzed in this explicit manner, even 

though the individual calculations to solve the Navier-Stokes equations or the structural dynamics 
equations are implicit. This procedure is described in detail, along with many examples, in 15,16 
The Navier-Stokes code should be reliable for fully attached and mildly separated flows. 


CLASSICAL INCOMPRESSIBLE POTENTIAL FLOW MODEL 


The aerodynamic loads input to Equation (14) was also calculated using the classical 
incompressible potential flow theory (linear aerodynamic theory) 1^. For an arbitrary time 
dependent motion, assuming the airfoil starts from rest at t = 0, 


Lift = 7 tpb 2 [h + Va - ba d] - 27tpVb 


wo.75c(0) d(s) + 


f 


dwo,75 C (q) 

da 


<J)(s-a) da 


(16) 


Moment = 7tpb 2 [bah - Vb(l -a) a - b 2 (1 + a 2 ) a 


- 27tpVb 2 (l- + a) 


wo.75c(0) <t>(s) 



dwp 75 c(q) 
da 


<Ks-a) da 


(17) 


where <()(s) is Wagner function, s = Vt/b and wo. 75c = -Va - h - b(^- a) a. Substituting h and t 
respectively for h and t and writing Cl, Cm in terms of Lift and Moment 



(18) 


C L = 7C 


h 

-2M 2 


+ 


a 

2M 


. 1 


( S 

a a 
4M 2 

- 2k 

w(0)<p(s)+ I ^-<(>(s-a)da 
Jo da 


c„=5 


ah 

L2M 2 


(i-a)d 

■%r 


(| + a 2 ) a 
4M 2 . 


2k (i + a) 


w(0) <j)(s) + j 4>(s-o) da 


(19) 


where w = wa 75 c/V = -a-^-A-(i--a)^ and s = ^ . It should be noted that in above 

equations, Mach number appears due to the choice made for nondimensionalizing time and 
velocity. It has nothing to do with compressibility effects and the flow is still incompressible. 
Equations (14), (18) and (19) are solved in the same fashion as Equations (14) and (15) for the 
same initial conditions to analyze the fluid-structure interaction. This aerodynamic model should be 
reliable for attached flows at small angles of attack. It is used to provide a reference for the other 
aerodynamic models. 


RESULTS AND DISCUSSION 

The flutter behavior of a NACA 0012 airfoil at a Mach number of 0.3 was studied with the 
proposed semi-empirical model, the Navier-Stokes solver and the classical incompressible potential 
flow theory. The particular flow conditions and structural parameters were selected, so that a 
comparison can be made to some of the earlier flutter studies 15,16 

For the flutter calculation using the semi-empirical model, first Equations (14) and (15) 
were solved with initial conditions as a = 15°, h = 0, Cl and as their steady state values at 15° 

angle of attack, and the first derivatives with respect to nondimensional time as 0. At these 
conditions, the airfoil is on the verge of stall. It was then released and was allowed to follow the 
motion dictated by the structural dynamics equations and the fluid dynamics model. The structural 
parameters used in this study are a^ = - 0.5, co^/g^ = 0.2, x a = 0.25, r a = 0.5, fi = 100, £h> 

= 0. Distances are nondimensionalized with respect to semi-chord (b), see Figure 3 for 
* 

definitions. The dimensionless speed V was varied between 4 and 8. The airfoil response to these 
two dimensionless speed values are shown in Figure 4. It was found that at V =4, the airfoil 
returned to steady state following a period of damped oscillations. For the higher flutter speed 

ifi 

case, V =8, the airfoil oscillations grew rapidly with time, and flutter was triggered, as clearly 



seen in Figure 4. As discussed subsequently, for some V between 4 and 8 the motion was a limit 
cycle flutter. 

The calculations shown above were repeated at lower angles of attack ( a between 0° and 
10°) and a similar flutter behavior was observed. The dimensionless speed V was varied between 
4 and 8. The airfoil response to these two dimensionless speed values at 5° angle of attack are 
shown in Figure 5. 

This was further analyzed by running the Navier-Stokes solver for these conditions and the 
same flutter behavior was observed. The starting point was the steady viscous flow over a NACA 
0012 airfoil at 0.3 Mach number and 4 million Reynolds number. At 15° angle of attack, the airfoil 
response are shown in Figure 6. Flutter can be seen at the higher speed and the damped oscillations 
at the lower speed. These calculations were repeated for 5° angle of attack and the same flutter 
behavior was observed. These are shown in Figure 7. 

The airfoil response were also calculated the classical potential flow model. The initial 
conditions for Equations (14), (18) and (19) were same as those for the semi-empirical model as 
explained earlier. At 5° angle of attack, the results are shown in Figure 8. It can be seen that flutter 
occurs at higher speed and damped oscillations at lower speed. These calculations were repeated 
for higher angles of attack and a similar flutter phenomenon was observed. Of course, using the 
classical potential flow (linear) model, the flutter condition should be independent of initial 
conditions and that is found to be so. 

For all three aerodynamic models, the damped and the flutter motions occurred at 
frequencies that roughly correspond to the natural frequency associated with the pitching motion. 
This suggests that the pitching mode is the most dominant mode in the flutter phenomenon 
considered here. In the present work, the Navier-Stokes solver was used to calculate the flutter 
response in fully attached flow (at a of 5°) and in mildly separated flow (at a of 15°). It was 
suggested in that a different type of airfoil response was observed in these cases, i.e. no flutter 
occurred in fully attached flow and flutter occurred in separated flow. In the present study, 
however, flutter occurred under both flow conditions for all three aerodynamic models. The semi- 
empirical model for prediction of aerodynamic forces works well for fully attached flows and 
mildly separated flows, but was not extended to massively separated flows as no data were 
available for fitting in that region. The flutter behavior shown by the classical potential flow model, 
which assumes an inviscid flow and considers no stall effects, agrees reasonably well with that 



predicted by the other two methods. The agreement among these three different procedures further 
confirms this flutter as essentially conventional flutter and not stall flutter, perse. 

Flutter boundary as a function of initial angle of attack is plotted in Figure 9. It can be 
observed that for angles less than 10°, a particular value of Vp* gives neutrally stable oscillations. 

For V* greater than this value, growing oscillations occur and for V* less than this value, decaying 
oscillations occur. The flow over airfoil is fully attached for these initial angles of attack. 

Note the velocity for onset of flutter decreases as angle of attack increases as intuition might 
suggest. However, at sufficiently large angle of attack, once flutter begins at some velocity it is a 
finite amplitude limit cycle oscillation until some larger velocity is exceeded when the oscillation 
then grows (exponentially) with time. For angles greater than 10°, neutrally stable, limit cycle 
oscillations occur over a range of values of V* instead of one particular value. Below this range, 
decaying oscillations are seen and above this range, growing oscillations occur. These oscillations 
are shown in Figure 10. A limit cycle behavior in the phase plane is observed over a range of V* as 
seen in Figure 1 1. The emergence of this limit cycle region coincides with the transition from fully 
attached flow to mildly separated flow over the airfoil. This behavior was also exhibited by the 
Navier-Stokes solver for similar initial conditions as shown in Figures 12 and 13. 

It should be noted, however, that the range of V over which the limit cycle is observed is 
somewhat different for the two fluid dynamics models. Typically the difference is modest, of the 
order of 20% in flow velocity. For the semi-empirical fluid model, a limit cycle of finite amplitude 
is observed over the range of V* ~ 5.25 - 6.75. For the Navier-Stokes fluid model, a limit cycle 

^ jjC 

was observed at V = 4.75, but not at V < 4 or V > 5.5. A comparison of results from the 
Navier-Stokes solver and the semi-empirical model for the maximum limit cycle amplitude is 
shown in Figure 14. The results are similar qualitatively although there are obvious quantitative 
differences. 

Because of the substantial computational cost, only a limited number of results were 
obtained using the Navier-Stokes solver. The flutter calculation using the Navier-Stokes solver 
took about 6000 seconds of computer time on a CRAY-XMP. Whereas, the other two methods 
required 4-5 seconds of computer time for a similar calculation. This shows the considerable 
savings in computer time (and memory), that can be achieved by use of the semi-empirical model 
in a standard aeroelastic analysis. Unfortunately there are no available experimental flutter data to 
compare to the present results from either the semi-empirical model or the Navier-Stokes model at 
high angles of attack. Note that the flutter results obtained from classical, linear potential theory are 
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included in the present discussion primarily as reference data. One would not expect the latter result 
to be meaningful at large angles of attack of course. However the classical aerodynamic theory has 
been experimentally validated for flutter at small angles of attack. 


CONCLUSIONS 


A system of ordinary differential equations can be used to model the aerodynamic forces on 
an airfoil for arbitrary motions. This semi-empirical model can be fitted for a variety of flow 
conditions including separated flows, provided experimental (physical or numerical) data exist for 
those conditions. This model then can be used with the structural dynamics equations to analyze 
flutter behavior of an airfoil in a standard aeroelastic analysis including conditions of mild flow 
separation. The savings in computer time and memory requirement are significant in comparison to 
the methods based directly on the Navier-Stokes equations. 

An interesting phenomenon is observed, where the structural-fluid interaction shows a limit 
cycle behavior as a function of V . For fully attached flows, this limit cycle behavior is seen for 
only a certain value of V*, whereas for mildly separated flows, it is seen over a range of V*. 
Further studies are being done to understand this phenomenon better. 

At present, this semi-empirical aerodynamic model is postulated, but efforts are underway 
to demonstrate its deduction from the first principles of fluid mechanics. 
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APPENDIX 1 

PARAMETERS FOR THE SEMI-EMPIRICAL MODEL 

(NACA 0012 airfoil at = 0.3 and Re = 4 million) 

Parameters for Equation (1) : 

AO / co A 2 = 0, Ai / co A 2 = 5.4878, A2 / co A 2 = 14.0997, A / 0) A 2 = -55.198. 

Ca = CaO + CaI a + Ca 2 « 2 ’ CaO = 2 - 898 > CaI = 54 - 701 > Ca 2 = -291.0861. 

“A = ^AO + ®A1 a + “A2 a2 ’ ®A0 = 7.0293, (0 A1 = -52.212, co A2 = 279.09. 
B A = b A 0 + B A1 a + b A 2 a2, b A0 = 242.709, B A1 = 415.509, B A2 = -2230.9. 

Parameters for Equation (10) : 

For a < 10.0174 \ 

A0 / co A 2 = 0, Ai / ® A 2 = -0.01557, A2 / 0) A 2 = -0.0998, A3 / co A 2 = 1.1753. 
For 10.0174° < a < 13.5636°, 

A0 / © A 2 = 1.466, Ai / co A 2 = -22.585, A2 / co A 2 = 1 15.66, A3 / co A 2 = -196.84. 
For 13.5636° < a , 

Ao / co A 2 = 4.314, Ai / co A 2 = -41.328, A2 / (0 A 2 = 129.49, A3 / to A 2 = -135.48. 
For all a, 

^A = CaO + CaI a » CaO = 3 - 4 004> CaI = 3 -3174. 

Wa = co AO + co A1 a , co A0 = 10.1120, o) A1 = -4.0498. 
b A = b A0 + B A1 « , B A0 = -71.87, B A1 = 27.9352. 
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Figure 1 . — Lift coefficient from semi-empirical model and 
Navier-Stokes solver (NACA 001 2 airfoil, a = 10° -5°cos(wt), 
M = 0.3, Re =4 million). 



Figure 2. — Moment coefficient from semi-empirical model and 
Navier-Stokes solver (NACA 001 2 airfoil, a = 1 0° -5° cos(wt), 
M = 0.3, Re = 4 million). 



Figure 3. — Two degrees-of-freedom typical section model. 
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Figure 8.— Time response of a two degree-of-freedom solid-fluid system using the classical incompressible potential flow theory 
(NACA 001 2 airfoil, M = 0.3, Re = 4 million, initial angle of attack = 5°). 



Initial angle of attack (deg) 

Figure 9.— Flutter boundary as a function of initial angle of 
attack. 
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(c) V* = 7.5. (f)V*=:7.5. 

Figure 1 0. — Time response of a two degree-of-freedom solid- Figure 1 0.— Concluded, 

fluid system using the semi-empirical model over a range cf 
V* (N ACA 001 2 airfoil, M = 0.3, Re = 4 million, initial angle of 
attack = 15°). 



(c) V* = 7.5. 

Figure 1 1 Limit cycle behavior in the phase plane for a two 
degree-of-freedom solid-fluid system using the semi-empirical 
model over a range of V* (NACA 001 2 airfoil, M = 0.3, Re = 4 
million, initial angle of attack = 15°). 


(f) V* = 7.5. 
Figure 1 1 .—Concluded 






a (rad) 


Figure 1 2. — Time response of a two degree -of- freedom solid-fluid 
system using the Navier-Stokes solver over a range of V* (NACA 
001 2 airfoil, M « 0.3, Re = 4 million, initial angle of attack = 1 5°). 


Figure 13. — Limit cycle behavior in the phase plane fora two 
degree-of-freedom solid-fluid system using the Navier-Stokes 
solver over a range of V* (NACA 001 2 airfoil, M = 0.3, Re = 4 
million, initial angle of attack = 15°). 


— Semi-empirical model 
0 Navier-Stokes solver 
A Linear Theory 


8 10 


Figure 14. — Maximum limit cycle amplitude as a function of V* 
(NACA 001 2 airfoil, M = 0.3, Re = 4 million, initial angle of 
attack =15°). 
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